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Development of a Zonal CFD Code for Hypersonic Flows 


INTRODUCTION: 

As detailed models of the National Aerospace Plane become avail- 
able, CFD codes will be required to accurately model flows over 
those models at Mach numbers up to 25. No current CFD code can 
handle both the high Mach number conditions and the complex 
geometries which will be required. 

Codes do exist which can produce Navier-Stokes solutions about 
simplified body shapes in the hypersonic regime. Four such codes 
have recently been compared by Edwards, et al., using a bicone 
shape and a generic elliptical all-body vehicled). While all of 
these codes are fairly accurate, they will become increasingly 
difficult to apply as model geometries become more complex. 

TNS (Transonic Navier-Stokes) is a CFD code which has been 
successfully applied to geometries as complex as the F-16A (2). 
The zonal capability of TNS allows a complex body to be modeled 
as an assemblage of relatively simple zones, with appropriate 
boundary conditions between zones. TNS was developed for 
transonic flow simulation over bodies which are best modeled 
using several zones, either due to the complexity of the geometry 
or to central memory limitations. The flow solver it uses is a 
diagonalized Beam-Warming method (3). While results in the 
transonic regime are very good, attempts to apply TNS to 
hypersonic flows have shown that this flow solver is probably a 
poor choice for hypersonic cases. 

This situation suggests that a synthesis of the TNS geometry 
capability with a more appropriate flow solver might produce a 
code capable of modeling detailed hypersonic vehicle designs. 

One flow solver selected for testing is that of F3D(4), which 
has been applied to STS at a Mach number of 7.9 and to a generic 
wing-body hypersonic configuration at Mach 25 (5). F3D was also 
included in the four-code comparison mentioned above. 

The work described here is based on a version of F3D provided by 
Y. Rizk, and a version of TNS which was modified for high-sweep 
cases by Jeff Cordova. 


TNS Testing: 

The existing TNS code has been tested to determine the highest 
Mach number at which solutions can be obtained. 

A generic elliptical all-body configuration was selected as a 
test case, since it represents the overall shape of a hypersonic 
vehicle in a simple way. A 2-D elliptic grid generator developed 
by Cordova was used to create Y-Z cross sectional grids about the 
body, and the grids were simply stacked in the X (streamwise) 
direction to create a three-dimensional grid. A more optimal 
grid could have been generated with a true three-dimensional 
method, but this was not done, since the flow solver will be 
required to tolerate some imperfections in grid orthogonality and 
smoothness as models become more complex. 


Once the grid was obtained, tests were made at increasing Mach 
numbers. Table 1 lists the cases which were run. The number of 
iterations shown reflects the time required to select effective 
time step and dissipation parameters, and may be much greater 
than the number of steps which would be required for another run 
of the same case. In all cases, the Reynolds number was 0.75E+06. 

For Mach numbers up to about 4, it proved fairly easy to obtain 
solutions. Each run was started from freestream conditions. A 
moderate increase in the fourth-order artificial dissipation was 
sufficient to maintain stability in the early steps, and the 
dissipation could be reduced to 0.01 or less after one order of 
magnitude convergence had been achieved. This value is the same 
as is used in transonic cases. 

When the Mach number reached 6, the dissipation parameter had to 
be set more carefully. It was found that the stability of the 
code peaked sharply when the correct dissipation was used (Figure 
1). Unfortunately, the correct value varied with Mach number and 
with the number of time steps which had been made, so the 
optimization must be repeated several times for each case. 

At a Mach number of 7.4, there was no longer a stability peak for 
any positive value of the dissipation parameter (Figure 1). The 
solution could no longer be obtained from free stream conditions. 
Instead, the Mach 6 solution was used as a starting solution, and 
the outer boundary conditions were ramped slowly up to Mach 7.4 . 
Solutions for positive angles of attack were obtained by similar 
slow changes in the boundary conditions. At 4 degrees, a 
solution was obtained, but at 6 degrees an instability in the 
wake region prevented convergence. 

Solutions were not obtained for Mach numbers greater than 7.4, 
even at zero angle of attack. While such solutions are not 
impossible, the increasing computational effort required seems 
prohibitive . 


CNS Development: 

In order to obtain solutions at higher Mach numbers, and to 
improve shock capturing in both supersonic and hypersonic 
conditions, a new code is being developed. The new code will be 
called CNS (Compressible Navier-Stokes ) . It will retain the 
zonal geometry capability of TNS, while incorporating a flow 
solver more suited to high Mach number cases. 

One of the flow solvers under consideration for CNS is the 
partially flux-split flow solver of F3D. The solver provides 
upwinding along one coordinate direction. In order to test this 
option, the F3D solver has been included as an option in the 
existing high-sweep TNS code. So far, no solutions have been 
generated with the new code. 


RESULTS : 


The solutions which have been obtained with TNS on the all-body 
do not correspond in Mach number, Reynolds number, and angle of 
attack to any published results obtained by other methods. 
Edwards, et al.(l) provide comparisons of results from four codes 
which have been validated on a bicone configuration and applied 
to the all-body geometry. To date this is the only data 
available, since no experimental results have been released for 
publication. Extensive validation will begin only after a flow 
solver has been incorporated which can give solutions over a 
wider range of Mach numbers. 

Figures 2 and 3 are included as examples of the data to be used 
for validation. Figure 2 shows boundary layer profiles at 3 X 
stations along the body centerline (symmetry plane) the 
conditions are M - 7.4, Re - 0.75e+06, alpha - 0 degrees. Figure 
3 shows the pressure coefficient along the centerline, above and 
below the body. The conditions are the same as for Figure 2, 
except at 4 degrees angle of attack. 


CONCLUSIONS: 

While the zonal grid system of TNS provides excellent 
modeling of complex geometries, improved shock capturing, and a 
higher Mach number range will be required if flows about 
hypersonic aircraft are to be modeled accurately. 

A new CFD code, CNS (Compressible Navier-Stokes ) , is under 
development to combine the required high Mach number capability 
with the existing TNS geometry capability. One of several 
candidate flow solvers for inclusion in CNS is that of F3D . This 
upwinding flow solver promises improved shock capturing, and more 
accurate hypersonic solutions solutions overall, compared to the 
solver currently used in TNS. 
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Table 1 


Cases tested using the Pulliam-Chaussee diagonalized algorithm in TNS 
with the all-body configuration. 


Freestream 

Alpha 

Grid Size 

Iterations 

Convergence 

Mach Number 

(degrees) 

( total ) 


(orders of 
in L2 norm 
inner zones 

magnitude drop 
of the residual) 
outer zones 

7.4 

6.0 

114,324 

735 

0.2-2. 2 

0.8-3. 6 

7.4 

4.0 

114,324 

1200 

3.0 

2.-5. 

7.4 

4.0 

114,324 

1702 

4.0 

3.-6. 

7.4 

2.0 

114,324 

900 

1.5-2. 6 

0.0-3. 9 

7.4 

0.0 

114,324 

3000 

4.-6. 

7. 

7.4 

0.0 

114,324 

115 

- 

- 

6.0 

0.0 

114,324 

2300 

4.-6. 

2.-5. 

4.0 

0.0 

348,595 

1200 

5. 

2. 4-5.2 

2.0 

2.0 

348,595 

1550 

4.5 

2. 3-4. 4 

2.0 

2.0 

114,324 

900 

3. 9-4. 5 

2. -3. 6 



Stability Optimization 




Boundary Layer Profiles 
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FIGURE 2 

Values at 3 streamwise stations on the centerline 
Mach 7.4, alpha = 0 degrees, Re = 0.75E+06 



Pressure Coefficients 



FIGURE 3 

Values along the all— body centerline 
Mach 7.4, alvha = 4 dearees, Re = 0.75E+06 





ORIGINAL PAGE IS 
OF POOR QUALITY 


Figure 5 


Standard Euler Set. numerical b.c.. Round-edged Wing. 64X64 Cell 
M^-2.0, a*10 u , P“70 g , e 2 »0.25, e 4 *0.01, 1. Surface pressure, 2. 

Crossflow Mach number, 3. Crossflow Velocity. 




Figure g. Zero-Total -Pressure-Loss Euler. Set ( 1 ) K 
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Figure 7. Standard Euler Set, C.f. b.c., round-edged wing, 64X64 cell, 
M w *2.0, a»10°, 6*70°, e 2 «0.25, e.»0.01, 1. Surface pressure, 2 
Crossflow Mach number, 3. Crossflow velocity. 
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Figure 8. Zero-Total -Pressure-Loss Euler. Set (1), C-f. b.c., round-edged 
wing. 64X64 cell, M.,-2.0. o*10°, p-70 u . e,«0.25, e,=0.01, 1. 

Surface pressure, 2. Crossflow Mach number, 3. Cfossflow Velocity 
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Figure 9. Standard Euler Set, numerical b.c., round-edged wing, 96X128 cell 
M^-2.0. a»10 u , P*70°, e 2 *0.25, e^-0.01,1. Surface pressure, 2. 

Crossflow Mach number, 3. Crossflow velocity. 
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Figure 10. Zero-Total -Pressure-Loss Euler, Set U). 

wing, 96X128 cell, M w «2.0. ot-10°, 3-70°, 
pressure, 2, Crossflow Mach number, 3, 


numerical b.c.. round-edged 
e 2 =0.25, e 4 =0.01, 1. Surface 

Crossflow velocity. 





Figure 11. Standard Euler Set, C.f. b.c., round-edged wing, 96X128 cell, 

M .*2.0, owlO 0 , 3*70°. e 2 *0.25, e^sO.Ol, 1 . Surface pressure, 2. Cross 
flow Mach number, 3. Crossflow velocity. 





Figure 12. Zero-Total -Pressure-Loss Euler, Set (1). c. f. b.c., round-edged wing 
96X128 cell, M^.O, a-10°, 0*70°, e 2 «0.25, e.-O.Ol. 

LSurface pressure, 2. Crossflow Mach number, 3. Crossflow velocity. 


o ittie tm 
O lifted • 
0 irtte WB 
A Ilf tend 





mm 


Will 


Hi! lift 


original; page is 

(OF POOR QUALITY 


Figure 14. Three-Oleenslonal Transonic Flow. Standard Euler Set. Sharp-edged delta 
wing, 80X38X48 cell, H.-0.7, a*15°. AR-1.5. Cj-0.12. €,-0.005. 

(a.,h.) 1. Surface Pressure, 2. Static Pressure Contours. 3. Crossflow 

velocity . , 

(<•.*•)!. Static pressure Contours, 2- Crossflow velocity 
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Figure 14. Three-D1«ens1onel Transonic Flow, Standard Euler Sot, Sharp-edged delta 
wing. 80X38X48 cell, M.-0.7, *.15°. ARM. 5. c 2 -0.12, e^-0.005, 

(•*8.) 1. Surface Pressure. 2 , Static Pressure Contours, 3. Crossflow 

velocity 

(«•,*)!. static pressure Contours, 2. Crossflow velocity 



Figure 15. Three-Dimensional Transonic Flow, Standard Euler Set, Sharp-edged delta 
wing. 80X38X48 cell. M,«0.7 t «-15°. AR-1.5, e 2 «0.12, c 4 «0.005. 
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Figure 16. Three-dimensional subsonic flow, Isentroplc Euler set, sharp-edged 


delta wing, 80x38x48 cell, M. • 0.3, a « 20.5°, AR « 1, c 2 ■ 
0.12, c H • 0.005, (a.,b.) 1. surface pressure, 2. static presure 
contours, Cn, 3. crossflow velocity 








Numerical Examples ( 2 ) 


flat- plate sharp-edged delta wing 

• Symmetric Conical FlmfCUssicJ e 0 ! cr e ^ s .) 

M. * 2, a * 10° and 0 * 70° 

• Uniformly Rolling Wing a Conical Flow (cl«ss;c-f E-Ur o>*-1 

•* 

w a Jf- “ 0.5 , M w » 2, a * 0° and 0 » 70° 

• Rolling Oscillation of a Wing In a Locally- Conical Flow (cU*».uF e -U r E r ».) 

« ■ “ « 0 cos kt e x , and 0 « - 9 (nax sin kt 
where 

0 max * w o /k ' K “ or * 

u -0.35,9 , n/12. k * 1-337, M. - 2 , <x -10 c and 0 - 70° 

0 " ax -t =■ 4.699 
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Fig. 1 Steady symmetric flow around a delta 
wing, M_=2, a=10°, 0=70° 



2. cross-flow velocity 



Fig. 1 Steady symmetric flow around a delta 
wing, H.-2. «»10® , 8-70° 
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Fig. 2. Uniform rolling of » delta wing. 
9 M./2. a*0°. 0*70°. u»0.$ 



4. static pressure 
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Fig. 3 Roll angle, angular speed and angular 

acceleration of the rolling oscillation 
motion. 
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Fig. 4 .Rolling oscillation of a delta wing. 
M -"2. 0-70°. u=0.35, k*1.337, 

W 15 °* t-0-1.07. 9=0-(-1455°) 
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2. cross-flow velocity 


1. surface pressure 


Fig. 4 . Rolling oscillation of a dalta wing. 
H.*2. a*10°. 0-70° , u-0.35. k»l.337, 

9 max* I5 °> 8*0-(-l<fl5°) 




3. cross-flow Mach 


4. static pressure 
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2 . cross-flow velocity 
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Fig. S. Rolling oscillation of a delta wing. 
M_-2. a*10°, 0-70°, w-0.35, k-1.337, 
8* *15° , t«l.07*-2.16. 8-H4.85 0 )- 
1-3.76°) 



3. cross-flow Mach 


4. static pressure 




ig. 6. Rolling osculation of a delta wing. 

«*10°, 0=70° , w=0.35, k-1.337. 
ta2 -16 -3.19. 0«(-3.76°)- 
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cross-flow Mach 


4. static pressure 




3. cross-flow Mach 


4. static pressure 











Fig. 9. Rolling oscillation of a delta wino 
M *2. a=10° fl» 7 n 0 a wing, 

. . « 10 , P*70 . w*0.35, k*1.337, 

0 nax xl5 .» t*5.35 -6.46, © = (-11.46° )- (-10.63°) 



1. surface pressure 
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Fig. 9 . Rolling oscillation of a delta wing, 
M„»2. «*10° , e*70°, u-0.35, k-1.337. 
9 *15°. t»S.3S*-W6, 9.(-U.46«)- 

(^0.63°) 



3. cross-flow Mach 4 - static pressure 






3. cross-flow Macrt 


4. static pressure 



















FIGURE 14- TIME HISTORY OF THE LIFT AND ROLLING-MOMENT COEFFICIENTS 
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CONCLUDING REMARKS 


1. ME HAVE PRESENTED APPLICATIONS OF THE CLASSICAL AND ZERO-TOTAL-PRESSURE-LOSS SETS 
OF EULER EQUATIONS TO SHARP- AND ROUND-EDGED DELTA WINGS- 

2- WE HAVE EXPLAINED THE ORIGIN OF THE TOTAL PRESSURE LOSS IN THE CLASSICAL SET . 

3- FOR SHARP-EDGED DELTA WINGS, ALL SETS OF EULER EQUATIONS PRODUCE THE SAME 
SEPARATED FLOW SOLUTIONS- THE SOLUTIONS REPRESENT REAL FLOW SOLUTIONS- 

4- FOR ROUND-EDGED DELTA WIN6S AND FOR COARSE GRIDS, THE SOLUTION DEPLNDS ON (1) THE 
LEVEL OF DISSIPATION, (2) THE ACCURACY OF THE SURFACE BOUNDARY CONDITION, AND (3) 
THE TYPE OF EULER EOS- SET- ONE CAN GET SEPARATED, PARTIALLY SEPARATED, OR 
ATTACHED FLOW SOLUTIONS- 

5- FOR ROUND-EDGED DELTA WINGS AND FOR FINE GRIDS, ATTACHED FLOW SOLUTIONS ARE 
OBTAINED. 

6- WE HAVE ALSO PRESENTED THREE-DIMENSIONAL FLOW SOLUTIONS AND ASYMMETRIC FLOW 
SOLUTIONS INCLUDING UNSTEADY FLOWS FOR SHARP-EDGED DELTA WINGS- 

7- EULER EOS- SHOULD BE RESTRICTED TO SHARP-EDGED WINGS FOR REAL FLOW SOLUTIONS- 
FOR ROUND-EDGED WINGS, NAVIER-STOKES EOS- MUST BE USED- 


